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Abstract 

We consider solutions of the Aw-Rascle model for traffic flow fulfill- 
ing a constraint on the flux at a; = 0. Two different kinds of solutions 
are proposed: at a; = the first one conserves both the number of 
vehicles and the generalized momentum, while the second one con- 
serves only the number of cars. We study the invariant domains for 
these solutions and we compare the two Riemann solvers in terms of 
total variation of relevant quantities. Finally we construct ad hoc finite 
volume numerical schemes to compute these solutions. 

Key Words: Aw-Rascle model, traffic models, unilateral constraint, Rie- 
mann problem, finite volume numerical scheme. 
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1 Introduction 

The paper deals with solutions to the Aw-Rascle vehicular traffic model [2] 

f dtp + d^{pv) = 0, . . 

I dty + {yv) = 0, ^ ' > 

satisfying a constraint on the first component of the flux at x = 0: 

p{t,f))v{t,Q)<q, (1.2) 

where g > is a given constant. Here p, v and y denote respectively the 
density, the average speed and a generalized momentum of cars in a road. 
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Moreover y = p{v+p{p)), where p S C^([0, +cx3[; [0, +oo[) is a pressure 
function satisfying 

p(0) = 0, 

p'{p) > for every p> 0, (1.3) 
p 1-7' pp{p) is strictly convex. 

Problem (jl.ip . ()1.2p models the presence of a constraint on traffic flow at 
the point x = 0, such as a toll gate, a traffic light, a construction site, 
etc. All these situations limit the flow at a specific location along the road. 
Conservation laws with unilateral constraints as ()1.2p were first introduced 
in [7], see also [HEIIP] for further analytical results and applications. In these 
papers, the scalar Lighthill-Whitham [17J and Richards |21j traffic model is 
coupled with a (possibly time-dependent) constraint on the flow, as in ()1.2p . 

The model presented here constitutes the first example of a system of 
two equations with constrained flux. The Aw-Rascle model (jl.ip belongs 
to the so-called "second order" traffic models, i.e. models consisting in two 
equations (see [U [201 [l2j for other examples). System (jl.ip can also be 
written 

dtp + d^{pv) = 0, 

dt{p{v+p{p))) + d^{pv{v+p{p))) =0. ^ • ^ 

The first equation in ()1.4p states the conservation of the number of vehicles, 
moving with flow rate pv. The second equation is derived from the former 
one and from the evolution equation of the quantity w = v + p{p) (often 
referred to as "Lagrangian marker"), which moves with velocity v: 

dt{v+p{p)) + vdx{v + p{p)) = 0. 



The system in conservative form (|1.4p belongs to the Temple class |22j . 
i.e. systems for which shock and rarefaction curves in the unknowns' space 
coincide. In particular, for such systems the interaction of two waves of the 
same family can only give rise to a wave of the same family. 

The Aw-Rascle model (|1.4p has been widely studied in the mathematical 
literature. Concerning the model itself, various extensions have been pro- 
posed, see O m [121 [131 [E]. The model can also be used to describe traffic 
flow on a road network, as explained in [TTl [T^ [T5]. 

In this paper we restrict the analysis to the Riemann problem for (jl.ip . 
()1.2p . i.e. to the Cauchy problem with piecewise constant initial data of the 
form 

{p\y^), if X < 0, 
{p\y^), ifx>0. 

We propose two Riemann solvers, described in Sections 2.1 and 2.2: the 
first one conserves at a; = both the number of cars and the generalized 
momentum, while the second one does not conserve the generalized momen- 
tum. In particular, the first Riemann solver produces a non-entropic shock 



(p,y)(o,x) 
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wave at X = 0, which travels with zero velocity. In Section 3 we describe the 
invariant domains corresponding to the two Riemann solvers, and in Section 
4 we compare the total variation of relevant quantities. Section 5 is devoted 
to the construction of ad hoc numerical schemes designed to capture the 
proposed solutions. 



2 The Riemann problem 

In this section we deal with the Riemann problem 

dtp + d^{pv) = 0, 

dt{p{v + pip))) + d^{pv{v + p{p))) = 0, 

/ N/n ^ / (P''^''), if 2; < 0, 
^^'")(0'") = i {p^,v% ifx>0, 



(2.1) 



in the domain T) = x , and with the constraint ()1.2p . 

We denote by f{p, v) the flux for system (|1.4p . and with fi{p, v), /2(/5, v) 
its components, i.e. 

For reader's comfort, we resume in the following tables the relevant quanti- 
ties concerning systems (jl.ip . (|1.4p respectively. In {p,y) plane they write: 

Ai = -p{p) + ^ - pp'ip) A2 = -p{p) + ^ 



I + pp'ip) 



n = I _J j r2 = 

VAi • n = 2p'ip) + pp"ip) > VA2 • ra = 

Liip;po,yo) = ^p ^2(^:^0,^0) = ^P + pipip) - Pipo)) 

z = -p- pip) w = ^ 

In (p, v) plane their expression is: 

Xi=v- pp'ip) \2=V 

-I \ f I 



''=^p'ip)) ''~-\^ 
VAi • ri = 2p'ip) + pp"ip) > VA2 • r2 = 

Liip;po,VQ) = VQ+pipo) - pip) L2ip;po,vo) = vq 

z = V w = V + pip) 

Above, Ai and A2 denote the eigenvalues of the Jacobian matrix Df, ri and 
r2 the corresponding right eigenvectors, Li and L2 the first and the second 
Lax curve, z and w the 1- and 2-Riemann invariant respectively. 

We remark that the system is strictly hyperbolic away from p = (i.e. 
Ai < A2). Moreover the first characteristic speed is genuinely nonlinear, 
with characteristic speed that can change sign, and the second one is linearly 
degenerate with strictly positive speed. 
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Definition 2.1 A Riemann solver for system i2. 1\) is a function, which 
associates, for every initial condition {p\v^) G V, {p'^,v^) ^ V, a map 
belonging to L-^(M) and representing a solution to 112. 1]) at time t = 1. 

By TZS we denote the classical Riemann solver for (|2.ip . i.e. the Riemann 
solver without the constraint ()1.2p : see for example [2]. We introduce some 
more notation. 

Given {p\v^) € V and g > 0, let us consider the set 

h = {pG [0,+oo[ : pLi(p;p',i;') =g} (2.3) 
= {p G [0, +oo[ : p(y' + p(p') - p{p)) = g} . 

The set Ii contains the densities of all the points (p, f) G P belonging 
to the Lax curve of the first family passing through (p' , ) and such that 
fi{p\v^) = q. If Ii 7^ 0, then we denote by p, v, pi, vi respectively 

p = max/i, v = ^, pi = min/i, vi = —; (2-4) 
P Pi 

see Figured) Given {p^ ,v'^) G T) and (7 > 0, let p2 and V2 be defined by 

P2L2{p2;p'',v'') = q, V2 = (2.5) 

P2 

i.e. {p2,V2) belongs to the Lax curve of the second family passing through 
{p'^,v^) and satisfies fi{p2iV2) = In particular, note that V2 = and 
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Given {p\v^) and {p'^,v^) £ V, let us consider the set 



h = {pe[0,+(X)[ : Li{p;p^,v^) = L2{p;p'',v'')^ 
= {pe[0,+^[:v^+p{p^)-p{p)=v'^ 



(2.6) 



and define 



max/21 v"^ = 



(2.7) 



which provide the intermediate state for the classical solution to ()2.ip . 
Lemma 2.1 Let (/>', u'), (p*", 7;^) G V and q > be fixed. Assume lll.3\) 



then the set Ii is not empty and it consists in exactly two different points: 
h = {/5i,p}. 

Proof. Notice that the function p 1— )• pLi{p; p\v^) is strictly concave by 
the hypotheses (|1.3p on the pressure function p{p) and so, by (j2.3|) . the car- 
dinality of 1 1 is at most 2. 

Denote with {p^ ,v^) the trace of 7^5((p^ ?;'),(/, ^^'")) at x = 0+. Since 
the waves of the second family have strictly positive speed, then v^^ = 
Li{p^ ] p\v^). Therefore, if Ii = or it contains only one element, then 
pM^M < ^ therefore fi{TZS{{p\v^),{p'' ,v~')){<d)) < q, which IS a con- 
tradiction. Thus the only possibility is that Ii is composed exactly by two 
elements. □ 

We propose two different ways of solving problem (|2.1|) - (|1.2|) . 

2.1 The Constrained Riemann Solver TZSl 

In this part, we introduce the Riemann solver 7^5^ for (j2.1|) - ()1.2p . which is 
characterized by the conservation of both the quantities p and y = p{v+p{p)) 
at X = 0. 

Fix {p\v^), {p^,v^) S P. The Riemann solver 7^5^ is defined as follows. 



1. If fi{nS {{p\v'), {p'\v-)) (0)) < g, then 

ns\ ({p\v'),{p\v^)) {x) = ns({p\v'),{p^y)) (x) (2.8) 



holds. If 



fi{ns{{p',v'),{p^,vnm)>q, 



for every x G M. 



2. If /l(7^5((p^^;'),(p^^;0) (0)) > q, then 



nSl[ip\v^),{p^,v^)) (x) = | 




if X < 0, 
if X > 0. 



(2.9) 



5 



Proposition 2.1 The Riemann solver TZS\ satisfies 
for every (/9^u') and {p^,v^). 

The proof follows directly from the construction of the Riemann solver T^iS^. 

Remark 1 The Riemann solver TZS\ is determined by imposing the con- 
servation of both quantities p and y at x = and by respecting the constraint 
condition see also for an example of a Riemann solver at a node, 

which conserves both p and y. 

In the following, we denote by w{TiS\ {{p\v^), (p^ v^) i^)) the w com- 
ponent of the Riemann solver TZSl (^{p\v^), {p^^v"^)) (x). 

Proposition 2.2 The Riemann solver TZSl satisfies the maximum principle 
on the second Riemann invariant w = v +p{p), i-c. 



mm 



< w{TlS\ (^{p\v''),{p^ ,v^)^ (x)) < max Mx 



The property easy follows from the maximum principle satisfied by the 
classical Riemann solvers IZS ((p', w'), (p, v)) for x < and IZS ((pi, ui), {p^ ,v^)) 
for X > 0, and by the fact that w = wi = wK 



2.2 The constrained Riemann solver TZS^ 

In this part we describe the Riemann solver 7^52' which conserves only the 
car density p at x = 0. 

Fix {p\v^) G [p'\v^') € T>. The Riemann solver TZS^ is defined as 
follows. 

1. If /l(7^5((p^^;0,(p^fO) (0)) < g, then we put 

nsl {{p\v%{p\v^)) {x) = ns({p\v'),{p''y)) (x) (2.10) 

for every x € M. 

2. If /l(7^5((p^^;'),(p^^;0) (0)) > q. then 

nsl ({p\v%{p^y)) (X) = j ^t^r '''^;?V*^l5u ^ '5 " ; n 

J IZS {{p2,v2)Ap' ,v''))[x), if X > 0. 



Proposition 2.3 The Riemann solver TZS2 satisfies 

fi{nsl{ip\v^),{p\v'-))m<q 

for every {p\v^) and {p^,v^). 



(2.11) 
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The proof follows directly from the construction of the Riemann solver 
TZSl 

Remark 2 The Riemann solver TZS\ conserves only the density at x = 0; 
therefore it is in the same spirit of Riemann solvers introduced for traffic at 
junctions in fllf - 



3 Invariant domains for TZSl ^'^l 

In this section, we want to describe the invariant regions for the Aw-Rascle 
system with constraints. First, we recall that, for every < vi < V2, < 
wi < W2 and V2 < W2, the set 

T^vi,v2,wi,w2 = {{P, v) e V : vi < V < V2, wi < V + p{p) < W2} (3.1) 

is invariant for (jl.ip : see Figure [2] and [16]. The hypothesis V2 < W2 implies 
that the Riemann invariants w = W2 and z = f 2 intersect in P at a point 
different from the origin. For a given g > 0, we define the function of class 

C2(]0,+CX)[) 

/i, :]0,+oo[ R 

V > — > p{^) +v, 

which gives the value of the Riemann invariant w of the point (p, u) € P 
such that pv = q. Indeed we have that hg{v) = w if and only if w = v +p{p) 
with pv = q. 

Lemma 3.1 Fix q > and assume There exists v = v{q) > such 

that the function hq{v) is strictly decreasing in ]0,v[ and strictly increasing 
in ]v, +oo\. 



Proof. We have 



h''iv) = ^ hp' ( + ( 1 



V^ L ' \V/ V 



and so, by (jl.3p . we deduce that /ig(i') > for every v > 0; this means that 
h'g{v) is a strictly increasing function. Note also that ()1.3p implies that 

lim p{p) = +00. (3-3) 

p— >+oo 

Indeed, if (|3.3p does not hold, then there exists M > such that p{p) < M 
and so pp{p) < Mp for every p > 0. This is not possible since the map 
p I—)- pp{p) is strictly convex. This implies that 

lim haiv) = +00 and lim h'(v) = 1; 

hence, since h'^ is a strictly increasing function, there exists a unique v > Q 
such that h'q{v) = 0. Therefore hq is strictly decreasing in ]0,'u[ and strictly 
increasing in ]v, +oo[. □ 
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Fi gure 2: The invariant domain 1^vi,v2,wi,w2- 

Proposition 3.1 Fix < vi < V2, < wi < W2, V2 < W2 and q > 0. If 

hq{v) > W2 for every v £ [^1,^2]; then I'j;i,t,2,w)i,«)2 invariant for both the 
Riemann solvers TZS\ and TZS\. 

Proof. The liypothesis hq{v) > W2 for every v € [^1,^2] implies that 
sup{/i(/0,w) : {p,v) E 'D^^^y2,wi,w2,v S [vi,V2],v + p{p) = W2} < q 
and so 

sup{/i(p,v) : {p,v) € Vy^^^2^u)i,w2} < q- 

Therefore the Riemann solvers TZSl and 7^52 in the domain 'Dy-^^jj2,wi,w2 co- 
incide with TZS. □ 



3.1 The Riemann solver TZSl 

The next proposition describes the invariant domains for 7^5f . 

Proposition 3.2 Fix < vi < V2, < wi < W2, V2 < W2 and q > 0. 

Assume and that there exists v G [fi,f2] such that hq{v) < W2- The 

set 2?^i,i,2,«)i,i«2 ^-^ invariant for the Riemann solver TZSl, if and only if 

hq{vi) > W2 and hq{v2) > W2; (3.4) 

see Figure\^ 



Proof. Clearly, if condition (|3.4p holds, then the set I^i,i,„2,«)i,j«2 is invariant 
for TZSl, since both {p, v) and {pi,vi) belong to 2?,)i,j;2,«)i,«)2 every possible 



choice of initial conditions in I'i,i,t,2,«)i,w2- 
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Assume now that I'«i,j;2,i«i,«>2 invariant for TZSl- 
If hq{vi) < W2, then denote with {p\v^) = {p^ ,v^) € I'j;i,,)2,«'i,«'2 the solution 
to the system 

f +p{p^) = W2, 
\ = Vl. 

By hypotheses, we deduce that /o'u' > q and so the trace of the Riemann 
solver TZS\{{p\v^),{p^ ^v"^)) at the point 0— is given by {p,v), which does 
not belong to I'i,i,t,2,«)i,w2! since hg{vi) < W2- This argument shows that 

hq{vi) > W2. 

If hq{v2) < W2, then denote with (p', v'') = {p^ ,v''") € 'Dy-^^v2,uji,w2 the solution 
to the system 

f +p{p^) = W2, 
\ v'- = V2- 

By hypotheses, we deduce that p'l;' > q and so the trace of the Riemann 
solver TZS\{{p\v^), {p^ ,v^)) at the point 0+ is given by {pi,vi), which does 
not belong to I'i,i,t,2,u;i,«)2> since hq{v2) < W2- This argument shows that 
hq{v2) > W2- This completes the proof. □ 



3.2 The Riemann solver USl 

In this part, we describe the invariant domains for T^cSg. First let us intro- 
duce the following necessary conditions. 

Lemma 3.2 Fix < < W2, < wi < W2, V2 < W2 and q > 0. As- 
sume lil.3\) and that there exists v £ [vi,V2] such that hq{v) < W2- If the set 
T^vi,v2,wi,w2 is invariant for the Riemann solver 71.82, ^hen hg{vi) > W2- 
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p 

Figure 4: The invariant domain I'«i,i;2,«)i,«)2 ^^^^ Riemann solver TZS\. 

Proof. Assume by contradiction that hq{vi) < W2- Denote {p\v^) = 
€ Vy^^v2,wi,w2 the solution to the system 

f +p{d) = W2, 
\ = V. 

By hypotheses, we deduce that p^v^ > q and so the trace of the Riemann 
solver TZS^iip'' ,v^), {p^ ,v^)) at the point 0— is given by {p,v), which does 
not belong to ^?i,i,„2,«i'i,w2' since hg(vi) < W2- □ 



Lemma 3.3 Fix d < vi < V2, ^ < wi < W2, V2 < W2 and q > 0. As- 
sume il.3\) and that there exists v G [^1,^2] such that hq{v) < W2- If the set 
T^vi,v2,wi,w2 ^-5 invariant for the Riemann solver TZS\, then hq{v) > wi for 
every v G [^1,^2]- 

Proof. Assume by contradiction that hq{v) < wi for some v G [fi,'L'2]- 
Denote {p\v^) = {p'',v^) G I^i;i,t)2,«>i,«>2 the solution to the system 

J +p{p^) = U)2, 
\ v'- = V. 

By hypotheses, we deduce that p^v'' > q and so the trace of the Riemann 
solver TZS'^Hp^v''), (p^jV^)) at the point 0+ is given by {p2,V2), which does 
not belong to I^j;i,,,2,?i>i,i«2 5 since V2 = v and hq{v) < wi. □ 

We have the following proposition about necessary and sufficient condi- 
tions for a domain to be invariant for 7^52. 
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Proposition 3.3 Fix < vi < V2, < wi < W2, V2 < W2 and g > 0. 
Assume \1.3\) and that there exists v € [vi,V2] such that hg{v) < W2- The 
set ^?t)i,?;2,i«i,u;2 invariant for the Riemann solver TZS^ (see Figure^ if 
and only if 

hqiyi) > W2 and hq{v) >wi G [wi,W2]- (3-5) 

Proof. By Lemma 13.21 and Lemma 13. 3^ we need to prove that condi- 
tion ()3.5p is sufficient in order I'^i,i;2,u;i,«)2 be invariant for the Riemann 
solver 7^52. Thus we assume that condition (j3.5p holds. 

Since 2^1,1 ,i)2,wi,w2 is invariant for (jl.ip . it is sufficient to prove that the left 
and the right traces at x = for 7^52 belong to I'„i,,)2,«'i,«)2- So fix {p\v^) 
and {p^ ,v^) in 'D^^^y^^wx,w2- '^S\{{p\v^),{p'^ ^v"^)) produces the classical 
solution, then we conclude. Assume therefore that lZS\{{p\v^),{p'^ ,v^)) 
does not produce the classical solution and denote with {p,v) and {p2,V2) 
the left and right traces at x = for TZS\{{p\v^), {p^ ,v^)). 
If {p,v) 'T>vi,v2,wi,w2i then every point {p,v) of the Lax curve of the ffist 
family through {p\v^) contained in P«i,t;2,«)i,««2 has the property that pv < 
q and so the Riemann solver gives the classical solution, since waves of 
the second family have strictly positive speed. This permits to prove that 
,■0) € T-^vi,v2,'wi,w2- 

If {p2,V2) ^ 'T>vi,v2,wi,w2^ then every point {p,v) of the Lax curve of the 
second family through (p'^,v^) contained in I^i;i,„2,wi,w2 has the property 
that pv < q and so the Riemann solver gives the classical solution. In fact, 
if > v^ , then a shock wave of the ffist family with strictly negative speed 
appears, if v^ = v^, then no wave of the first family appears, whereas if 
v'' < , then all the states {p,v) of the rarefaction wave have flux pv less 
than or equal to q. This permits to prove that {p2,V2) € 'T'vi,v2,wi,w2- 

The proof is thus completed. □ 



4 Total variation estimates for IZSf and TZSl 

In this section we make a comparison between the two Riemann solvers TZSl 
and 7^52 in terms of the changes in the total variation of various quantities. 

Fix (p', u'), (p'', w*") € v. We denote with pi and p2 respectively the p- 
component of TZS\{{p\v^), {p^,v''')) and of TZS\{{p\v^), {p'^,v^)). Moreover 
we denote with vi, V2 respectively the w-component of TZS\{{p\v^), {p^jv"^)) 
and of TZSl{{p'-,v'-),{p'',v'')). Finally, we put yi = pi{vi + p{pi)), m = 

P2{V2 +P{P2))i = ^'l +PiPl), W2 = V2 +p{p2)- 

In order to facilitate the reading of the following calculation, we refer to 
Figures and El 
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pv 



L2 



pv_ = q- 



P™P2 Pi 



P"P2 Pi P P 

Figure 5: Notations used in the paper: case pi > p2- 



y 





Pi P2 P"" P P pi P2 p^ p p 

Figure 6: Notations used in the paper: case p\ < p2- 



4.1 Total variation of the density p 

This subsection deals with Tot.Var.(pi) and Tot.Var.(p2)- The foUowing 
proposition holds. 

Proposition 4.1 For every initial conditions {p\v''), {p^ ,v^) G V, we have 
that 

Tot.Var.{pi) > Tot.Var.{p2). (4.1) 

Proof. linS\{{p\v^),{p'\v')) = 7^5^((p^^;'), (p^-u'')), then Tot.Var.(pi) = 
Tot.Var.(p2)- Therefore, we assume that 



nsw,v^l{p\v^))^nsl{{p\v%{p\v^)). 

In this case we have that fi{TlS{{p\v''), {p'^ ,v^)){0)) > q and so, by con- 
struction of TZS^ and TZS^ , we deduce that pi{x) = p2{x) for a.e. x < 0. 

Moreover, for x > 0, Tot.Var. ^P2|]o+cx,[) ~ \ f>2~ p'^\i since the states {p2,V2) 
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and {p^ , v'^) are connected by a contact discontinuity wave of the second fam- 
ily. Hence 



Tot.Var.(p2) = Tot.Var. P2| ) + \p-p2\ + \h - P 



P 



+ \P- P2\ + \P2 



First consider the case = Li{p'^ ; p\v^), so that {pi,vi) and {p^,v^) can 
be connected by a wave of the first family. We get that 

Tot.Var.(pi) = Tot.Var. + |p - pi| + |pi - /| 

+ \p- pi\ + \pi- p1- 

If p'^ < Pi, then p2 G]p^,Pi] and we obtain Tot.Var. (pi) = Tot.Var. (p2). 
If p"^ > pi, then pi < P2 < p^ < p and so Tot.Var. (pi) — Tot.Var. (p2) = 
2(P2 - pi) > 0. 

Consider now the case v'^ ^ Li{p^; p\v^). We have that 



Tot.Var. (pi 



p -p 



+ ip-pii + ipi-p'"i + ip--pn. 



If P2 < Pi, then we get that p"* < p2 and 

Tot.Var.(pi) - Tot.Var.(p2) = p2 - p"^ + |p"^ - p'\ - \p2 - pH > 
by the triangular inequality. If P2 > Pi, then we get that p™ > p2 and 



Tot.Var.(pi) - Tot.Var.(p2) 



P2+p'"-2pi + ip--pn 



\P2 



> 2(p2-pi)>0 
by the triangular inequality. This completes the proof. 



□ 



4.2 Total variation of the velocity v (i.e. the first Riemann 
invariant) 

This subsection deals with the total variation of the velocity, i.e. of the first 
Riemann invariant z. 

Proposition 4.2 For every initial conditions {p\v''), {p^ ,v^) G V, we have 
that 

Tot. Var. {vi ) > Tot. Var. {v2 ) • (4.2) 

Proof. If 7^<Sf((p^'^;')> (P^^^0) = T^Sl{ip\v^),ip'' ,v'')), then the thesis 
clearly holds. Therefore we assume that 

nsW,v'),{p\vn)7^nsli{p',v^),ip^,vn). 
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In this situation we have that /i(7?.5((p', w'), (p^, t;''))(0)) > q and so, by 
construction of 7lS\{{p\v^),{p'^ ,v'')) and 'RS\{{p\v^),(p^ ^v"^)), we deduce 
that vi[x) = V2[x) for a.e. x < 0. It is clear that 



Tot.Var.({;2) 



since V2 = L2{p2; p"^ ,v^) = v^- 
If v"^ = Li{p^; p\v^), then 

Tot.Var.(i)i) = 



V — V 



V — V 



+ \V - V2\ + \V2-V 



+ \v — V 



V — V 



+ \v — Vl \ + \vi — v'^ 



and so, by the triangular inequality, we deduce Tot.Var.(i;i) > Tot.Var.({;2). 

If / Li{p^; p\v^), then 



Tot.Var.(?)i) 



V — V 



V — V 



+ \v - Vl \ + \vi - V \ + \v 
+ \v — Vl\ + \vi — \ 



since f*" = v'^ by ()2.7p . Again, Tot.Var.('i;i) > Tot.Var.(i)2) by the triangular 
inequality. 

The proof is so finished. □ 



4.3 Total variation of the generalized momentum y 

This subsection deals with the total variation of the generalized momentum 
y = p{v + p{p)). 

Proposition 4.3 Assume that hypothesis il.3\) holds. For every initial con- 
ditions {p\v^), {p'^,v^) G V, we have that 

Tot.Var.{yi) > Tot.Var.{y2). (4.3) 

Proof. lillS\{{p\v^),{p'' ,v'')) = 7e5^((p',t;'), thenTot.Var.(yi) = 

Tot.Var.(y2)- Therefore we assume that 

nsl{{p\v'), {p^v^)) / nsl{{p\ v% (/, v^)). 

In this situation we have that /i(7^5((p', f'), (p^, f'"))(0)) > gandsoyi(x) = 
y2{x) for a.e. a; < 0. We define y = p{v + p{p)), y^ = p\v^ + p{p^)), 
y'^ = p''{v"~ +p{p'')), yi = pi{vi +p{pi)), m = P2{v2 +p{p2))- Note that 
y > max{yi,y2}- 
We have that 



Tot.Var.(y2) 



-y 



+ \y- M + \y2-y 
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Consider first the case = Li{p^] p\v^), which imphes that 



Tot.Var.(yi) = - y 



+ \y- yi\ + \yi - /I • 

and consequently 



If < pi, then, by (jl.Sp . we easily get that y'" < 2/2 < yi 
~ " 'li) = Tot.Var.(y2). 

then, by ()1.3p . y^ > y2 > yi and so 



ii p ^ pi, tnen, oy [\i.6\) . wt 
Tot.Var.(yi) = Tot.Var.(y2). 
If / >pi 
2(y2 - yi)- 

Consider now the case 7^ Li{p'''; p\v^). We have that 



Tot.Var.(yi)-Tot.Var.(y2) = 



Tot.Var.(yi) = y^ - y 



where y™' = p'^{v^ +p{p^))- 
If Pi < P2, then, by ()1.3p . we deduce that yi < 1/2 



■ Pi 

Tot.Var.(y2) = 

< 
< 



+ |y-yil + |yi-y™l + |y"'-yn, 

^ y2 < and so 



2/ -y 
.1 



y - y 

If pi > P2, then, by ()1.3p . we 
Tot.Var.(y2) = 



+ (y - y2) + |y2 - y''| 
+ (y-yi) + |y2-y"| + |y"-yn 
+ (y - yi) + |yi - y"l + |y" - yH = Tot.Var.(yi) 

deduce that yi > y2 > y*^ 

+ (y - y2) + |y2 - y'" 



■■"^ and so 



y - y 



< 



The proof is completed. 



y' - y + (y - y2) + |y2 - y""! + ly"" - yi 
y' - y| + (y - y'") + ly'" - yH = Tot.Var.(yi). 



□ 



4.4 Total variation of the second Riemann invariant w 

This subsection deals with the total variation of the second Riemann coor- 
dinate w = V + p{p)- 

Proposition 4.4 For every initial conditions (p', w'), f^) € D, we have 
that 

Tot.Var.{wi) < Tot.Var.{w2). (4.4) 

Proof. imS\{{p\v^),{p'' ,v'-)) = 7^cS^((p^^;0, (p'^wO), thenTot.Var.(5i) = 
Tot.Var.(z2). Therefore we assume that 



ns\{{p\v%{p\v^-))^nswy),{p\v^)). 
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In this situation we have that /i (7^5 ((/?', f'), (p**, t'''))(0)) > (7 and so t()i(x) = 
W2{x) for a.e. x < 0. We define w = v+p{p), = v^+p{p^), = v^+p{p'^), 
wi = vi +p{pi), W2 = V2 +p{p2)- Note that = w = wi. 
We have that 



Tot.Var.(zZ;2) 



v/ — w 



+ \ w — W2 \ + \W2 — W 



Consider first the case v^' = Li{p'^; p ,v ), which imphes that 



Tot.Var.fiDi) 



w — w 



Consider now the case ^ Li{p'^; p\v^). In this case we have that 



< Tot.Var.(it;2)- 



Tot.Var.(tZ;i) 



w — w 



+ \w'- 



w 



where w"^ = +p{p'^). Since = w, we conchide by the triangular 
inequality. 

The proof is so finished. □ 



5 Numerical schemes 

This section is devoted to the construction of finite volume numerical schemes 
to capture the solutions corresponding to 7^5^ and TZSl- 

Let Ax and At be two constant increments for space and time discretiza- 
tion. We then define the mesh interfaces Xj^i/2 = j'Ax (so that xx/2 = 
corresponds to the constraint location) and the cell centers Xj = (j — 1/2) Ax 
for j € Z, the intermediate times = nAt for n S N, and at each time 
t" we denote u" an approximate mean value of the solution of (jl.ip . (jl.2p 
on the interval Cj = [xj^i/2,Xj^i/2)i j ^ ^- In other words, a piecewise 
constant approximation x — > u(t",x) of the conserved variables u = {p,y) 
is given by 

u(r, x) = for all x G Cj, j £ Z, n G N. 
When n = 0, we set 

1 f^j+1/2 

Uj = — — / uo(x)(ix, for all j G Z, (5.1) 

where uq = {po,yo) G P is a given initial data (we will restrict the study to 
Riemann-type initial data). 

Given a sequence (u")jgg at time t", we concentrate now on the com- 
putation of an approximate solution at the next time level t""*"^. 

We will concentrate on Godunov scheme and show how to adapt it in 
order to match the constraint condition (jl.2p at x = 0. We recall that, as 
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pointed out in [5] , classical conservative schemes (like Godunov method) may 
generate important non-physical oscillations near contact discontinuities. 
For this reason we will restrict to Riemann data lying on the same second 
Riemann invariant, i.e. we take + p{p^) = +p{p^)- More general cases 
can be treated for example combining the techniques presented here with 
the Transport-Equilibrium scheme described in [5]. Note that, in any case, 
a contact discontinuity appears when applying the Riemann solver 7^52- 
For sake of completeness, we recall that classical Godunov scheme writes 

= < - - for all J G Z, (5.2) 

where the numerical fluxes are given by 

f;+i/2 = f "i+i) = fiT^i^l Ui+i)(0)) for ah j G Z, (5.3) 
and the usual CFL condition 

^max{|A,(u,")|, i = l,2}<l (5.4) 

holds. In the following sections we describe how to modify the definition of 
the numerical flux (j5.3p for j = 0. The simulations have been performed 
taking p(p) = p and Ax = 0.002. 

5.1 The Constrained Godunov scheme for TISX 

We follow the idea introduced in for the scalar case. We redefine the 
numerical flux at the interface Xi/2 = to take into account the imposed 
constraint p.2|) . We denote by f"j_,_i/2' ^2*j+i/2 ^'^^ components of the clas- 
sical Godunov flux: 

V ^2,j+l/2 J 

For j = 0, we replace it by i^^^ where 

f"l/2 = min{f[^,/2 ' l} ' 

fu _ fn ^Illl-^^^lfn JIlIlX ^^-^^ 

^2,1/2 - h,i/2 7^ - ™™ \ h,l/2 ' Q ( ■ 

1,1/2 L 1,1/2 J 

We stress that the above construction preserves conservation, in agreement 
with the conservative character of TZSl- 



Theorem 5.1 (Maximum principle) Under the CFL restriction ([57 
the finite volume numerical scheme defined by ()5.2p . (|5.3p and (|5.5p satisfies 
the maximum principle property 

inf {vf +P(P?)) < v^+Pip"^) < max {vf + pip^)) 
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for all j £ Z and all n € N, where = ( \ and v'j +p{p^) = ^ = ""^j • 

Proof. We observe first that the above maximum principle property on 
the second Riemann invariant is satisfied by the classical Godunov scheme 
()5.2p . ()5.3p (see for example [5| Remark 3.1 (ii)] for a detailed computation). 
Thus we only need to check what happens for j = 0, 1. 

If ^11/2 = f"i/2' then also ^21/2 ~ ^21/2 scheme reduces to the 

classical Godunov scheme. Therefore we assume that — Q < ^11/2 

and ^2 1/2 ~ ^2'i/2/^"i/2 ^ ■^2'i/2' ^^^^ case, recalling the construction of 
TZSl, it is easy to see that 

f[^,/2 = flK,u) = fl(ui,U?), 

f£,/2 = g^7e5«,u?)(0))=f2K,u) = f2(ui,O, 

where u and ui are, respectively, the left and right traces at x = of 
TZSKuq, u"). In fact, since TZS (uq , u) counts only waves of negative speed, 
we have that 

f «, u) = / {TZS (u^, u) (0)) = /(u) = ( J ) ■ 

On the other side, since TZS (ui, Uj^) counts only waves of positive speed, we 
have that 

f(ui,O = /(7e5(ui,O(0)) = /(iii)= f J ) . 

Hence the following bounds hold for j = 0, 1: 

inf l^r, ii), t&il < tt;" < max {wp, t£), , 
/=-!,. ..,2 ■' i=-l,...,2 

and we conclude observing that w = 'W\ = Wq. □ 

We have tested our method on Riemann data lying on the same 1- 
Riemann invariant, in order to avoid spurious oscillations due to the presence 
of contact discontinuities. More general data can be dealt with using the 
technique presented in [5j. Figures [71 [HI shows that the numerical solu- 
tions are in good agreement with exact solutions. In particular, our scheme 
perfectly captures the nonclassical shock at x = 0. 

5.2 The Constrained Godunov scheme for 7^52 

The Constrained Riemann Solver TZS2 is not globally conservative at the 
point X = (by definition, conservation holds only for the first equation 
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density at t=0.1 with rhoL=l . 5, rhoR=l . 5 



elocity at t=0.1 with vL=3, vR=3 



-0.1 0.0 0.1 0.2 0.3 0.4 0.5 



Figure 7: Test la : Solution of the constrained Riemann solver TZS\ with 
data = p'^' = 1.5^ = v'^ = 3 and q = 3: exact solution (dashed line), 
numerical approximation (continuous line). 



density at t=D.l with rhoL=4, rhoR=l . 5 



ith vL=Q .5, vR=3 



-0.1 0.0 0.1 0.2 



-0.5 -0.4 -0.3 -0.2 -0.1 



Figure 8: Test lb : Solution of the constrained Riemann solver TZSl with 
data = 4, p*" = 1.5, = 0.5, = 3 and g = 3: exact solution (dashed 
line), numerical approximation (continuous line). 
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in (jl.ip and therefore only car density p is conserved). As a consequence, we 
look for a non-conservative numerical scheme, i.e. we define two numerical 
fluxes 7^ ^1/2' such that 

= Uo-|^(fr/2"-f-i/2), (5.6) 

= <-^(f3-/,-f5+). (5.7) 

We set 

en,— cn,+ ■ ( cn 1 

fi;i/2 = fi;i/2 = ™in|fi,i/2>9j ' 

and 

£•71 

Jn,- _ nn 2,1/2 
^2,1/2 ~ n,l/2 fn 

1,1/2 

In order to capture the right trace at x = 0, we could envisage using a ghost 
cell type method (introduced in |10j . see also pS] and references therein 
for other applications), computing the ghost value u" corresponding to u", 
whose (p, v) components are given by 

yTl 

where f " = — —p{Pi)- This is obtained using the following flux 
Pi 

f2"v2 = '?K+p(9K)) ' 

whenever f"^/2 ^ 1- Unfortunately, due to the convexity assumption ()1.3p 
on the function p ^ pp{p), the velocity component is overestimated during 
the projection step of Godunov scheme in (|5.7p (see [5]). Therefore, the right 
trace cannot be captured properly: the velocity component is overestimated 
and the density is underestimated, see Figures El [TTJ In fact, at each time- 
step, we have v'^'^'^ = v""*"^ > and p^'^^ < p", where the inequality is 
strict generally speaking. 

In order to overcome this difficulty, we propose to simply keep the value 
of the velocity component fixed for j = 1, i.e. to replace the value ob- 
tained by ()5.7p with i2'i/2 ~ ^21/2 v^^^ = , and then updating the 
conservative component as 

y^+'=pr' K+p(pr')) , 

whenever f"]^^2 ^ 1- This allows to capture precisely the right trace of the 
discontinuity at a; = 0, as shown by numerical simulation in Figures [9l [TTl 
Only, a small amplitude oscillation traveling at speed v = is produced, 
see Figures [lOl [l2j 
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density at t=0.1 with rhoL=l . 5, rhoR=l . 5 



elocity at t=0.1 with vL=3, vR=3 



1.3 O.i 0.5 



Figure 9: Test 2a: Solution of the constrained Riemann solver T^iSg with 
data /o' = = 1.5, v'' = = 3 and q = 3: exact solution (dashed line), 
ghost cell method (dash-dotted line), our method (continuous line). The 
rectangles select the zoomed areas plotted in Figure [TOl 




Figure 10: Test 2a: Detailed view of a part of the computational domain 
of Figure M 
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density at t=D.l with rhoL=4 , rhoR=l . 5 



-0.1 0.0 0.1 



.3 0.4 0.5 



ith vL=0.5, vR=3 



Figure 11: Test 2b: Solution of the constrained Riemann solver T^iSg with 
data (} = A, = 1.5, = 0.5, u'' = 3 and g = 3: exact solution (dashed 
line), ghost cell method (dash-dotted line), our method (continuous line). 
The rectangles select the zoomed areas plotted in Figure [T2l 




Figure 12: Test 2b: Detailed view of a part of the computational domain 
of Figure [JTJ 
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density at t=0 . 1 with rhoL=4 , rhoR=l . 5 



4 . 0- 



3 . 5- 



3 . 0- 



2 . 5- 



2 . 0- 



1 . 5- 



1 . 0- 




. 5- 



. o-| 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 



velocity at t=0 . 1 with vL=0.5, vR=3 



4.5n 



4.0- 



3 . 5- 



3 . 0- 



2 . 5- 



2 . 0- 



1 . 5- 



1 . 0- 



0.5- 




. o-i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 



density at t=0 . 1 with rhoL=l . 5 , rhoR=l . 5 



3 . 5- 



3 . 0- 



2 . 5- 



2 . 0- 



1 . 5- 



1 . 0- 



. 5- 



. 0-| 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 



velocity at t=0 . 1 with vL=3, vR=3 

4.5- 
4.0- 
3 . 5- 
3 . 0-- 
2 . 5- 
2 . 0- 
1 . 5- 
1 . 0- 
0.5- 

. oH 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

-0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 




